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THE CONSISTENT SECOND-ORDER THEORY 
OF WAVE/STRUCTURE INTERACTION 



by 

C. J. Garrison 



ABSTRACT 

The consistent second-order theory of the interaction of regular gravity 
waves with a fixed object in water of finite depth is developed. The theory is 
carried out for the most general case of a body of arbitrary shape which may 
extend through the free-surface or be completely immersed. The incident wave 
evolves in the development as a second-order Stokes' wave. Boundary-value 
problems are established for both the first- and second-order velocity poten- 
tials and a numerical method based on the Green's function is outlined. 

INTRODUCTION 

The determination of forces exerted by gravity waves on large structures 
immersed in the sea has become of great practical interest in recent years. For 
example, in the design of bottom-mounted oil storage facilities or large ocean 
caissons, the wave-induced horizontal and up-lift forces and overturning moments 
are factors of primary importance. The effect of large amplitude waves in parti- 
cular is of importance in the determination of the permanence of an ocean struc- 
ture and, therefore, a higher-order theory appears to have significant practical 
value. For example, Apelt and Macknight (1976) found measured forces on a ocean 
caisson model in fairly large-amplitude shallow-water waves to be considerably 
in excess of calculations based on linear diffraction theory. 

With this application in mind then, the solution through the second-order 
in wave height is developed herein for regular wave interaction with a fixed 
object of arbitrary shape in water of finite depth. The theoretical development 
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is carried out in a mathematically consistent manner in that all terms through 
the second-order in wave height are included in the analysis, and, as such, 
represents the extension of the linear theory presented previously by Garrison 
and Seetharama Rao (1971). The problem is recognized as a regular perturbation 
problem in the small parameter (wave height/characteristic body dimension) and 
the incident wave appears within this framework as the second-order Stokes' 
wave. A numerical method for solving the resulting equations suitable for digital 
computer evaluation is outlined. 

REVIEW OF LITERATURE 

The solution of the linear wave/structure interaction problem is now fairly 
wel 1 -developed for infinite depth as well as for the finite depth case. A 
number of papers have appeared in the literature dealing with the oscillation of 
two-dimensional bodies on a free surface in water of infinite depth, a problem 
which is mathematically similar to the fixed-body/wave interaction problem. 
Examples include the work of Ursell (1949), Porter (1960), Vugts (1968), and 
Paulling and Richardson (1962). Dean and Ursell (1959) made both an experimental 
and theoretical study of small amplitude wave interaction with a fixed, semi- 
submerged circular cylinder in deep water. Yu and Ursell (1961) treated the 
problem of a semi -submerged circular cylinder oscillating vertically in water of 
finite depth. 

Somewhat less attention has been given to solving the corresponding three- 
dimensional problems. Apparently the first analysis of the linear wave inter- 
action with a fixed body was carried out by Havelock( 1940) for the case of a pile 
in water of infinite depth and MacCamy and Fuchs (1954) solved the same problem 
for water of finite depth. Havelock (1959) also evaluated the added mass and 
damping coefficients for a heaving semi -immersed sphere in deep water, and Wang 
(1966) extended this to include the finite depth case. Haskind's relations 
(Haskinds (1957)) as discussed by Newman (1962) may be used to determine the 
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linear or first-order forces acting on the fixed body from a knowledge of the 
damping coefficients for the same body oscillating in still water. Garrison (1974) 
has made this evaluation using Havelock's results for the hemisphere and the re- 
sults were in good agreement with his direct calculation of the heave force based 
on the Green's function approach. Kim (1965, 1966) has applied the Green's 
function approach to compute the excitation forces, added mass and damping co- 
efficients for a semi -immersed ellipsoid in deep water. 

Most work on wave/body hydrodynamics reported in the literature treats con- 
figurations and conditions which have primary application to the study of ship 
motion; only recently an interest in large, fixed offshore structures has develop- 
ed. Even though, probably the first potential flow (or diffraction) solution to 
a wave force problem having direct application to large ocean structures was ob- 
tained over twenty years ago by MacCamy and Fuchs (1954), namely, the small 
amplitude wave interaction with a fixed vertical circular cylinder (pile) in 
water of finite depth. This diffraction solution represents the only closed form 
solution of its type available in the literature and, therefore, is of primary 
importance to this general problem. Little work of this type having application 
to large ocean structures has appeared since that time until recently when 
Garrison and Seetharama Rao (1971) applied the Green's function approach to calcul- 
ate the first-order pressure distribution and resulting forces acting on a bottom- 
mounted hemisphere. Their experimental results for this configuration compared 
well with the theory. Garrison and Chow (1972) then extended this analysis to 
include fixed bodies of arbitrary shape and compared their theoretical results 
with experimental results corresponding to two different submerged oil storage 
vessel configurations. Milgram and Halkyard (1971) also have developed a linear 
theory and applied their method to certain axisymmetric bodies in deep water. 
Garrison (1974, 1977) has more recently given details of a practical method 
based on the Green's function to evaluate added mass and damping coefficients for 
floating bodies of arbitrary shape in water of finite depth, and at the same time 
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calculate the linear wave forces. 



It may be stated that, in general, linear theory is quite well developed 
and good agreement between theory and experiment has been obtained for wave forces 
corresponding to small amplitude waves in both the two- and three-dimensional 
cases. This is true as well for bodies oscillating in a free surface in still 
water. Experience has also shown that linear theory gives good wave force results 
when the body is fairly deeply submerged even if the amplitude of the incident 
wave is not small. However, if the object is surface piercing or not deeply sub- 
merged and the water is shallow, nonlinear effects caused by waves of finite amp- 
litude become pronounced. In such cases linear theory is inadequate, and a higher 
order wave/structure interaction theory becomes necessary. 

Nonlinear potential solutions to wave/structure interaction problems are 
much more limited than linear solutions. Olgilvie (1963) solved for the first- 
and second-order forces on a horizontal, submerged circular cylinder in water of 
infinite depth, the wave crests being parallel to the cylinder axis. However, 
his second-order forces included only the time-average part and, as a consequence, 
the first-order potential only was needed and obtained. This work, therefore, did 
not actually represent a complete second-order solution. However, a consistent 
second-order theory for a horizontal cylinder in still water has been developed 
by Lee(1968) and Potash(1970) . Garrison and Smith(1977) have treated the corres- 
ponding wave/fixed-body interaction problem. 

Recently there has been some interest in developing analyses for nonlinear 
wave interaction with three dimensional fixed bodies in water of finite depth. 
Chakrabarti (1972) made an attempt at the extension of MacCamy and Fuch's (1954) 
linear solution for the vertical circular cylinder (pile) to the fifth order. 
However, he treated the free-surface boundary condition improperly so that the 
results are mathematically inconsistent and, consequently, of little value for 
making judgements regarding the magnitude of nonlinear diffraction effects. More 
recently Raman, et.al. (1976, 1977) have also treated this problem 
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although the method of solution is extremely complex and there appears to be no 
simple means of checking the results. 

Isaacson (1977) has argued that the second-order solution for a body with 
a circular waterline curve does not exist because the free-surface and kinematic 
boundary condition on the cylinder are not consistent at the mean waterline on 
the cylinder. However, this inconsistency does not preclude the existance of a 
solution; it appears that such inconsistencies at the point where two different 
boundary conditions join is rather common. 

In this report the solution through the second-order in wave height is 
developed for bodies of arbitrary shape held fixed in regular waves in water of 
finite depth. The problem is formulated as a regular perturbation problem, and 
the incident wave is shown to represent a second-order Stoke 's wave. It is also 
shown that the second-order problem is the same type of boundary-val ue problem 
as the first-order problem, the primary difference being that the second-order 
free surface boundary condition is nonhomogeneous . The solution is formulated 
for both the first- and second-order potentials in terms of a Green's function, 
and the numerical results are presented for several configurations. 



FORMULATION OF THE PROBLEM 

The problem under consideration is depicted in Figure 1. A rigid object 
having a characteristic dimension a is immersed in water of depth h, and 
a train of regular waves propagates in the positive x-direction. The poten- 
tial solution is sought to the interaction of the incident wave with 
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y = 0 plane : mean water elevation 



FIGURE 1 DEFINITIONS 
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the fixed, rigid body through the second-order in the wave amplitude, 
or more precisely, in the ratio of the wave amplitude to characteristic 
dimension of the object, a\ 

The fluid is assumed to be irrotational so that the velocity 
potential, $, may be defined as 

f = V<$>Cx, g,2,i) (!) 

where q" denotes the fluid velocity vector. Assuming the fluid to be 
incompressible, it follows that the velocity potential must satisfy 

= O (2) 



within the fluid region. (The barred quantities denote dimensional 
quanti ties. ) 

On the bottom the kinematic condition specifying zero normal 
velocity is 

) ~ O ( 3 ) 

where h denotes the mean fluid depth. On the rigid, wetted surface of 

the body, specified by S(x,y,l) - o , the kinematic boundary 

condition describing the zero normal velocity condition is given by 

Y§> * VS = O (4) 

The elevation of the free surface above the y = 0 plane is denoted 
by 3 (x,z,t) and the kinematic boundary condition. 
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= o 



(5) 



; fa! 




must be satisfied on y = n. Moreover, on the free surface the pressure 
is set equal to zero, and, accordingly, Bernoulli's equation supplies the 
dynamic free surface boundary condition, 

j,* y l) f ,2, iff ?/M f*- <£/*, b^f]* S }<■*'£,*)= 9 ^ (6 ) 

where h denotes the Bernoulli constant. 

It is convenient to cast the boundary value problem thus far established 
in dimensionless form, and for this purpose introduce the following dimen- 
sionless parameters: 

x = x/a , y= y/ai j z=z/a , d=d/a ) h. =k/a. 

H= h/zq. y if - h*/a , cr*-a/g , t = z r± y {!) 

<p = cr <£> /qa 



where a = 2tt/T, T being the wave period. The dimensionless time is denoted 
by t, and h* denotes a dimensionless Bernoulli constant. As indicated in 
Equation (7), all the coordinates and length scales are made dimensionless 
with the characteristic dimension of the fixed object, a\ 

Using these parameters, Eqs. (2-6) which define the boundary value 
problem may be rewritten concisely in dimensionless form as: 



M ,*/* ) * o 


in the fluid 


(8) 


4> y (*,-h,2^) = O 




(9) 
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<^ X <y, £ / £ ) = ° 



on 



SCx.y.Z) = O 



00 ) 





Equations (11) and (12) are applicable outside the body on the free surface 
when the body is surface piercing. 

PERTURBATION PROCEDURE 

According to the methods of perturbation analysis, the potential 
function, free surface elevation, and Bernoulli constant may be expanded, 
respectively, in the small parameter, e, as 



where e relates to the wave amplitude in a yet undetermined manner. 

It may be noted that in Eqs. (11) and (12) the surface elevation, | , 
and consequently, the small parameter, e, does not appear explicitly. 
Accordingly, it is necessary to expand the particular functions involved 
so as to form functions which do contain e explicitly. For example, 
substituting Eq. (14) into (13) and expanding each term in a Taylor 




(13) 



+■ JCx. z, £) + 



(14) 



/f = h* + Oc 6 s ; 



(15) 
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series about e = 0 gives: 




Similar results can be obtained for the derivatives of 4> which appear in 
Eqs. (11 ) and (12). 

Substituting Eqs. (13-16) as well as expressions similar to Eq. (16) 
for the derivatives into the boundary- value problem given in Eqs. (8-12), 
the separate boundary-value problems for <t> j , <j> 2 , etc. are obtained. The 
first- and second-order problems may be isolated and are given, respectively, 



by: 



Fi rst- order: 



v i 4?( x ,y- %,-t) =o 



07 ) 



cb(x,-h, 

b 



(18) 



d) (X, y, 2 , t) = O on Scx,l J; z )=0 
1 *1 
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- y } cx,z,i) — o 
'y '* 



( 20 ) 



— O 

1 't 



( 21 ) 



Second-order: 



y ,*,0 ~ o 



( 22 ) 




( 23 ) 
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= o 



on 



=0 



(24) 




+ 4>Cx,o, Z.-6) j (x,tL,i) 

'i ' x 



* '£ 





Both 4 > i and <j > 2 are, in addition, subject to a suitable radiation condition 
which limits the scattering disturbance to outgoing waves. 

A great deal of similarity is evident between the first- and second- 
order problems. In fact, the only basic difference lies in the right 
hand side of the free surface boundary conditions; in contrast to the 
first-order free surface boundary conditions, the second-order boundary 
conditions are not homogeneous, the right hand sides of which being 
dependent on the first-order potential. 

Before proceeding to the solution of the boundary-value problems 
developed, advantage may be taken of the fact that the first-order poten- 
tial will be represented by a periodic function. Accordingly, we may 
define the complex potential u(x,y,z) as 



where R e denotes the real part, and a = 2ira/L, a being the characteristic 
dimension of the body and L the wave length. The symbol b denotes an 
unknown real constant. It is also appropriate and common practice in 
linear interaction problems to express the potential 44 as the sum. 




(27) 
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<t> = 4> x + 4? 



(28) 



where <j> r denotes the incident wave potential, and 4^ denotes the scattering 
potential which is due to the presence of the rigid body. The complex 
potentials u* and u^ are then defined in relation to ^and 4>* according to 
the form of Eq. (27) so that it follows that 

U, = U* 1- U, S ( 2 g) 

It is recognized, moreover, that the incident wave potential, 4^, must 
satisfy the first-order problem when no body is present, i.e., 4^ must 
satisfy Eqs. (17), (18), (20), and (21). These equations represent simply 
the boundary-val ue problem for the well-known, first-order progressive 
wave. The solution to this in terms of the present notation is given by: 

af = L c oshfaCh± y) J e mx (30a) 

a c °sh(ak) 

In Eq. (30a) a is defined in terms of v by the familiar expression from 
linear wave theory 

V^CT z a/g = a fankiah) (30b) 

where h denotes the dimensionless water depth. 

Now, substituting Eqs. (30a), (29), and (27) into the first-order 
problem defined by Eqs. (17-21), and eliminating between Eqs. (20) and 
(21), the following boundary-value problem is established for the first- 
order complex scattering potential: 
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V z U^Cx,4,z) = o 



(31) 



u.! y (x-M) = o 

U^Cx,4,z) = - osh l ( . QLh) [n a sinhLP-Ch+y)] + 1 n x co shLaCh^y)]] e ax 

C>n 5tx, y, 2 ) - O 

U*(x, O, i) - y U*(x, 0,2) - O 

The components of the unit vector directed normal to the immersed 
surface into the fluid are defined as n - l n x + J r n l . i n 
addition to Eqs. (31-34), u^(x,y,z) is also subject to the usual 
radiation condition which allows only outgoing scattered waves at a 
great distance from the body. 

Before dealing with the solution to the first-order scattering problem, 
the second-order problem will first be established. Proceeding in this 
direction it is appropriate in view of the linearity to express the solu- 
tion to the second-order problem, Eqs. (22-26), as the sum 

\ = i * € (35) 

where, similar to the first order problem, <t> r 2 denotes the second-order 
incident wave potential and ^ denotes the scattering potential. 

We again envoke the condition that when no body is present, there is 
no scattered wave so that <j>^ = ^ = 0. Moreover, under these conditions 
the boundary conditions on the immersed surface, Eqs. (19) and (24), are 
not applicable. Accordingly, we may substitute ^ ^ and <p 2 = <t> r , along 

with the Eq. (30) and (27), into Eqs. (22), (23), (25), and (26) and 



(32) 

(33) 

(34) 
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obtain, upon eliminating n 2 between Eqs. (25) and (26), the following 
boundary-value problem for the second-order incident wave potential: 

V z (p\x,Lj,z,i) = O (36) 

cb 1 (x , -h H, £ ) = O (37) 

z y 

4>(x,o,z ,t) + V (p T (x,o,z ,-6 ) - - \ b 2 (a l - v l ) s'm[z(ax-t )] (38) 

The periodic solution to the boundary-value problem specified by 
Eqs. (36-38) may be expressed as 

= jL a b z y R e [ i- U*(x,Lj) <e ] ( 39 ) 

where the second-order complex potential, u*, is given by 

u[cy,tj) = - -L- <e c£ax ~ (40) 

** a 

This expression is familiar in wave theory and is of exactly the same 
form as that given by Eagleson and Dean (1966) for second-order Stokes' 
waves. 

To complete the establishment of the boundary-value problem for the 
second-order scattering potential, it is necessary to substitute the 
forms 4 > l = 4 ^ + <j>^ and <j > 2 = <j > 2 + <j>* into Eqs. (22-26) and eliminate n 2 
between Eqs. (25) and (26). In addition, the known expressions for the 
incident wave potentials, Eqs. (30) with Eqs. (29) and (27) for the first- 
order potentials, and for the second-order potentials, Eqs. (40) with 
Eqs. (39) and (35), are utilized to obtain: 
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v z 4(o<, y,* ; £ ; = o 



(41) 



4? Cx,-h,z,£) = o 

* U 



(42) 




In Eq. (44), U(x,z)'is a time independent function which is generated by 
substitution of <f> into Eqs. (25) and (26). However, since the function 
is not needed in the second-order theory for purposes of evaluating the 
pressures and forces, it is for brevity's sake not written out in full 
here. 

The boundary-value problem now established for as given in Eqs. 

(41-44) is linear, and except for U(x,z) occurring in Eq. (44), is time 
w i 2 1 

dependent like e . Accordingly, it is appropriate to express the 
solution in the form 



where the second term denotes the part of the complex potential which 
is independent of time. By substitution of Eq. (45) into the boundary- 
value problem described by Eq. (41-44), separate boundary-value problems 



<p* = 3 Q b z y R e [i- e LZt +■ i U s z (x,<j,z)] 



(45) 



S ^ 5 $ 

arise for u 2 and u 2 • However, as will be evident subsequently, ^ 



t 
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only is needed in order to evaluate the hydrodynamic pressures and 
resulting forces to the second order and, consequently, the time 
independent part of <j>!| is of no interest. He, therefore, dismiss it 
from further consideration and concentrate on the solution u^. 

Substituting Eq. (45) into Eq. (41-44) the following boundary 
value problem is established for the second-order scattering potential: 



= o (46) 

, z )=0 ( 47 ) 

- ~ SlA ^ ah) [ n d Sinh[ 2 .Q-CK+-y)] + i. IT* coshLzaCh+ y)]^ 0 L ' 2aX (48) 
UlUo,?-) - 4-y u*(x,o,2) = f(\z) ( 49 ) 



where 



■fu.z) = | Jl/^ uf s U? n * u'(u: n -SVU^) + ^ u* ( (50) 

< < - 3 <L 0 



- 4 u. 



The radiation condition limits the scattered waves to outgoing 
regular waves and is given by 

li** 1 -\/7i ( -l. - lq, ) // $ - r> 

“V ( 51 ) 

where a^ denotes the dimensionless wave number for the second-order 
problem which is defined by 

A- V* ~ (X v -fanll (q l h) (52) 
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INCIDENT HAVE 



At this juncture it is appropriate to completely define the incident 



(26) for nj and n 2> respectively, and substituting the results into the 
expression for the elevation of the free surface defined by Eq. (14) 
gi ves 



Equation (53) expresses the free surface elevation in terms of the total 

potentials and, consequently, includes the effect of the scattered wave 

as well as the incident wave. However, it is presently of interest to 

obtain an expression for the free surface elevation of the incident wave 

alone in the absence of the object. For this purpose, therefore, we set 
s s 

<p = <t> = 0 and evaluate Eq. (53 j using the known expressions for the 

1 2 

incident wave potentials, Eqs. (30) with (27) and Eq. (40) with (39). 

The result is: 



wave and specify the unknown constants b and h*. Solving Eqs. (21) and 




(53) 



fo.z.t) = b T&[ C 



t(<xx- i 




b z (y l -Q l ) . l* 

4-y 1 



] j H- 0<C>) 



(54) 



* Ott co'JUah) Ufc-oshizah) R» c 



i ?-(ax » *L ) 




birirt (ah) 
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where n r denotes the dimensionless surface elevation due to the inci- 
dent wave in the absence of the body. However, the constant term in 
Eq. (54) must vanish so 



u* _ bUoLz 

* 4-y 



( 55 ) 



The remaining terms in Eq. (54) represent sinusoidal variations, the 
second-order contribution having twice the frequency of the first. 

If we define the wave height as the elevation difference between 
the trough and crest of the incident wave, then eb must represent the 
wave amplitude since the second-order contribution to the surface 
elevation in Eq. (54 ) is the same at both the crest and trough, the 
difference being zero. Therefore, in terms of the dimensionless wave 
height as defined in Eq. (7) we have: 

€ b = H = H /Zb- (56) 



where H denotes the elevation difference between the crest and trough. 
The incident wave profile is then specified by the dimensionless 
surface elevation as: 

3c H cos(ax-t) + cosh (ah) + coshCzak)) (57) 

* sink (ah) 

. Cos[Z Cax-O] + OCH 3 ) 

Equation (57) agrees with the expression for the second-order Stokes' 
wave given by, for example, Eagleson and Dean (1966). 
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PRESSURE, FORCES AND MOMENT 

The forces and moments acting on the immersed surface are deter- 
mined by carrying out surface integrals of the pressure. For this 
purpose the pressure is obtained from Bernoulli's equation as 

^ />/<£_**- 4>~ + ] -f>sS + PS h * ( 53 ) 

This expression may be cast in dimensionless form by use of Eq. (7), 

(13), (14), (15), (27), (35), (45), and (56), and the result carried 
to the second-order in wave height. The resulting expression is then 

f (59) 




) e-‘ a J +\a„\\ 



U, 



K 




} 



The dimensionless pressure coefficient expressed by Eq. (59) is 
defined as 



r - 



p_ 

P 3 & 



(60) 



in which P denotes the fluid pressure. 

The first term in Eq. (59) denotes the hydrostatic pressure as 
y represents the dimensionless depth beneath the mean free surface. 

The second and third terms are harmonic, the second-order part having 
twice the frequency of the first-order part. The final second-order 
term in Eq. (59) is independent of time and gives rise to time-average 
or steady state components of force and moments. 

The components of the dimensionless force and the moment vectors 
may now be expressed in terms of integrals of the pressure over the 
wetted surface area. That is, we may write 
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QW = ff f q. ds 



Z 



= h 2, 



(61 ) 



where the dimensionless force coefficients associated with the 
x, y, and z directions are defined, respectively, as 



C ,(0 = 

Pd °- 3 



(62a) 



C^t) - —Ed±l (62b) 

= - n ^ C - 3 - (62c) 



in which F x (t), Fy(t), and F 2 (t) denote the three components of force. 
Similarly, the moment coefficients are defined as 




MM 

Pg Q* 



(62d) 



C s (l) = — (PkiAL 

pg^ 









( 62 e ) 



(62f) 



in which M x (t), My(t), and M z (t) denote the three components of the 
wave-induced moment. 
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The functions g-j are defined as 



3, = R * » 9^ = n, # g, - 



(63a) 



g 4 = Cd+ y) n z - zn y J g 5 =z^-xri 0 g & = x n y - Cy+d) n* (63b) 



The forces and moments in Eq. (62) are defined relative to a coordinate 
system parallel to 0(x,y,z) but shifted downward a distance d as is 
evident from Eq. (63b). The unit vector directed normal to the wetted 
surface is defined by rt= l + j + k r? 2 and dS denotes an 
element of surface area made dimensionless with the characteristic 
dimension of the body squared, a 2 . 

We may make an initial ordering of the force or moment coefficient 
as given in Eq. (61) in terms of the small parameter, e, or equivalent 
H, as 



and the coefficients, p 0 . Pi, and p 2 , are defined in an obvious manner 
by comparison of Eq. (65) with Eq. (59). 

It may be noted at this point that, if the body is completely 
submerged, the surface integrals in Eq. (64 ) refer simply to the actual 
surface area of the body which is, of course, specified. However, if 




where the pressure coefficient is expressed as 




(65) 
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the object is surface piercing, S denotes the actual wetted area which 
is dependent on the instantaneous surface elevation of the water on the 
body. Accordingly, it is necessary to express the limits of integration 

on the wetted surface S in terms of n and carry out the expansion retain- 

2 

ing terms with coefficients up through H . 

For this purpose the dimensionless differential surface area may be 
written as dS = dc dl where dc denotes a 

dimensionless differential arc length on the wetted surface in the 
horizontal plane, and dl denotes a dimensionless differential length 
in an orthogonal direction. Thus,. it is possible to 

express dl in terms of y so that 

d S — dccJi = ~ 4== ^ y : (65) 

V / - 

The integrals indicated in Eq. (64 ) are now carried out over the wetted 
surface with y running from y = -e to y = n. Accordingly, a typical 
integral in Eq. (64) ma y be expressed in the form 




t i + 




j± 

b l 



h' 



t°j dc dy 



y~/- /?/ 



y= -e 



(67) 



where e denotes the lower limit of the immersed surface. 

The integrals of the type given in Eq. (67) are ne xt expanded in a 
Taylor series about H = 0 using Leibniz's formula for purposes of 
evaluating the derivatives of the integrals. Carrying out these 

expansions, substituting the results into Eq. (64), and retaining terms 

2 

up to order H yields: 
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( 68 ) 




■h 




- 3kLh ^ g. kdc f (7(» 3 ) 

Z > V I ~ fT-y’' / y=o b y 



where S 0 denotes the surface area below the plane y = 0 and C 0 denotes 
the closed waterline curve formed by the intersection of the y = 0 
plane with the surface of the rigid object. 

The first term in Eq. (68) represents the hydrostatic force or 
moment on the body where the associated displaced volume is defined as 
that beneath the y = 0 plane. However, since our interest is in the 
dynamic forces and moments, we may disregard this first term with the 
understanding that the buoyant force defined as such should be accounted 
for in order to determine the final total force or moment. 

Finally, applying Eq. (21) and (27) as well as the definitions 
for p 0 , p L and p 2 we obtain: 



C L U) - Ha 



[f R'iu. e ‘ *1 9i ds +H ! jjgff F?4«q‘u L - u,: 






)e“*J 9 i ds 
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Equation (69>a) may also be written in coefficient form as 




t Cf(H 3 ) 



(69 b) 



where the first- and second-order force (or moment) coefficients and 
phase shift angles are defined by comparison of Eq. £9 a) with $9 b) 
as follows: 



The steady state (nonperiodic) force (or moment) coefficient is also defined 
by comparison of Eqs. (69 a) and (69b) as 



The first term in Eq. f>9 b) represents the linear solution having 
the fundamental frequency, a. The next term which represents a second- 
order contribution to the force or moment is also harmonic but at twice 
the fundamental frequency. The last second-order term represents the 
steady state contribution which is independent of time. This latter 




(69c) 




(69d) 




'O 



where the dimensionless force coefficients, Fj. and F 2 ., are real. 




2V /- nj 



I ) 9 . & 



(69 e) 
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force is generally referred to as a drift force in ship hydrodynamics. 
(See, for example, Maruo (I960).) However, drift force is generally 
calculated by applying the momentum equation to the diffracted wave 
rather than by using near field results. 

The evaluation of the force or moment coefficients defined by Eq. 
(33c-e) represents the primary object of this paper. However, in 
order to carry out the evaluations indicated, it is clear that, given 
the geometry of the rigid immersed object, it is necessary to evaluate 
Uj and u 2 as well as the derivatives of ui on the surface area denoted 
by S Q and along the waterline curve C Q . 

SOLUTION TO THE FIRST-ORDER PROBLEM 

Having now established the need for the potentials, u ] and u 2 , we 
return to the consideration of the solution to the boundary- value prob- 
lem developed. The practical method of solution of the first-order 
problem given in Eqs. (31-34) is discussed in depth by Garrison(1974,1978) 
and, therefore, only the major steps will be outlined here. Following 
the Green's function method of solution we write uf,-'as the integral 
over the immersed surface S„ as: 



where U,n,c) denotes points on the immersed surface, fj U,n,0 denotes 
the unknown source strength and dS = dS/a 2 denotes the differential 
surface area on the immersed surface made dimensionless with the charac- 
teristic dimension of the object, a. The function, G (x,y,z;£»n»4.\>), 
denotes the Green's function which must satisfy the equation 




o 



(70) 
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V 2 G(x,y,E; g,}, 5; V ) = 8<*-$> S( •j-f) &(*-*) 



(71 ) 



in which 6 denotes the Derac delta function. The function, G(x,y,z;? ,n,t*,v) 
must also satisfy Eqs. (32), (34), and the radiation condition. Such a 
function is given by Wehausen and Lai tone (1960) as 



G(x, y , 2 ; €, j, j; v > = +■ 



i 

R 



I 

R' 



? pv f Cu±3dl LoskRchtj^)] j- , 

* ^ PV J smh («k) - V Cosh^h) ^ d * 



(72a) 



/ 2 77 YO,*- vll coshjACj+hjl cosh[Qj(J±h)l j (q i/l -) 

<?h - >>*/? + ^ 



in which 

R =[(*-£?- + cy-}f+ (z-~sy-]' /l ~ 

R = [(x-if-hCy+z-h + if-,- (z-'s) 2 -]''* 

XL + "$ )*] ^ 

The symbol, aj, is defined in terms of h and v as the solution 



(72b) 

(72c) 

(72d) 
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of the equation 



Ci, -tanh (d 1 h'> - V - o 



(72 e) 



In view of Eqs. (72e) and (30b) the symbol a^'s clearly equivalent 
to a. However, in the case of the second-order Green's function, the 
equivalence does not hold and, therefore, separate notation is maintained. 



An alternate series form of the Green's function is also given by 
Wehausen and Lai tone (1960) as 



where J Q and Y 0 denote, respectively, Bessel functions of the first and 
second kind of order zero, and K Q denotes the modified Bessel function 
of the second kind of order zero. In Eq. (73a) aj is defined by Eq. 
(72e) and the cuanti tiesx^ are defined as real positive roots of the 
equation 



In Eq. (72a), P. V. denotes principal value of the integral. 





(73b) 



The solution to the boundary-value problem stated in Eqs. (31 -34) 
is given by Eq.(7^); it remains, however', to determine the source strength 



function in order to evaluate the potential. This is accomplished by 
application of the boundary condition, Eq. (33), which results in the 
following integral equation from which fj may be determined: 



47T 




-,IJ , z; S, h d. s 



.. <? 



l ax 



CoshCah) 



[ r?y s > nhLo .( h + y )]+ in * cosh [ a ( fui)]J (74 ) 



The normal derivative of G evaluated on the immersed- surface, as required 
in Eq. (74), is determined in a straightforward manner by differentiation 
of either Eq. (72a) or (73a). 



SOLUTION TO SECOND ORDER PROBLEM 

The solution to the second-order problem may be formulated by application 
of Green's theorem to the fluid region denoted by Y"in Figure 2. Applying 




FIGURE 2 Region for application of Green's theorem. 

Green's theorem with subjects u^ and (the second-order Green's function) 
and assuming that G 2 satisfies the equation, 

V l G z (X t y t Z;S,1,S) = &(x-$>8(ij-})&(d-S) (75) 

as well as the kinematic condition on S^ given by Eq . (47) , the radiation 
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condition given by Eq.(51) on and the homogeneous free surface con- 

dition, 

J,s) - ;},;,]!> - o (76) 

gives the result: 

u£(y.y,t)-~ <Z„<*'V*;UT)- <£ l /, j)K <V, J)J^ s 

s ° (77) 

- ~ v] c/s 

The first term in Eq . ( 77 ) may be considered to represent a source and 
doublet distribution over the immersed surface of the body. However, 

Lamb (1932) has shown for the case of an infinite fluid that the potential 
may be represented as either a source distribution only, a doublet distri- 
bution only, or some combination of both. This proof is easily extended to 
the present case allowing the first term to be written as a distribution of 
sources only. Then, if the free surface boundary condition, Eq.(49), is 
used in the second integral the result is 

u£er.y,t) = fjS,lT) c/s 

(73) 

- j-f ds 

S s 

in which is the unknown source strength distribution. 

The second-order scatter potential defined by Eq . ( 78 ) satisfies Eqs. 

(46, 47, 49-51). The remaining kinematic boundary condition on the surface 
of the body results in the integral equation 
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( 79 ) 







s 



which is to be satisfied for x,y,z on S c . 

In both Eq.(73) and (79) a practical interpretation of the limits of 
the free surface area, , should be understood. In deriving Eq.(78) 
through application of the Green's theorem the outer limit of S 5 was 
considered to be large enough that the radiation condition was satisfied 
on S,*, . in practical application this occurs rather rapidly. 

In view of the conditions which were placed on the Green's function in 
the above development is clearly very similar to the Green's function 
used in formulating the first-order problem. In fact, if V is replaced 
by 4V and is replaced by in either Eq . ( 72 ) or (73) the re- 

sulting function is the appropriate Green's function for solution of the 
second-order problem. That is. 



where G is given by Eq . ( 72 ) and (73) and tf^is defined through Eq.(72e) 



In the case of the alternate series form of 0>2 given in Eq.(73), the roots 
are defined through Eq.(73b) as 




(80) 



a z innh Q L h - 4V = 



(81a) 



4 V — O 



(81b) 
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An alternate form of Eq.(72a) for the sources on free surfaces is given by 





x± CTcC^a.) 



yU smh Mk - 4V cosh m/i 



(82) 



, i M cu cqshjg^jJ^ 

4- V h -t- 5iK)h^a z h 



This if the form of the Green's function which is obtained by setting |-o 
in Eq.(72a) and which was also given by Wehausen and Laitone (1960) for 
the case of a harmonic pressure distribution on the free surface. The 
nonhomogeneous form of the free surface boundary condition given in Eq.(49) 
is equivalent to the boundary condition for the first-order problem for a 
pressure distribution of the free surface. 

NUMERICAL PROCEDURE 

Eq.(69) clearly indicates that the determination of the forces and 

moment acting on the immersed object rests on the determination of the 

potentials u-j , and derivatives of u^ at points on the immersed surface. 

The first-order scattering potential is specified by Eq.(70) in terms of the 

first-order source strength function, f^ . Thus, to determine u-| , it is 

first necessary to solve the integral equation, Eq.(74), for f-| . 

We may develop a numerical solution to Eq.(74) beginning with the 

partitioning of the immersed surface S Q into N subdivisions or facets of 

area A S j , each with a nodal point at its center located at the point 

(x., y., z.). Recognizing, moreover, that f, (£,?,J) is a well-behaved 
J J J 

function for smooth bodies we may define 
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jf . = £ ( yj , y,- , Zj) 



(85) 



and accordingly we may approximate Eq . ( 74 ) by the complex matrix equation: 

(X,:(y) £; = ( 86 ) 

Once . . (V.) has been evaluated by use of Eq . ( 83 ) the solution of Eq . (86 ) 

^ J 

may be carried out to determine f j ^ at each nodal point on the surface of 
the body. A discussion of the details of the numerical evaluation of /0\ 7 , ; 
including the evaluation of the l/R- singularity, has been given by 
Garrison (1978). 

Having determined the source strength function, the first-order poten- 
tial and its derivatives on the surface of the object may be determined from 
Eq . ( 70 ) . For this purpose we may replace the surface integrals with 
summations, writing 



< = p.x»> £. , 

"4 = i 




l;J = i, i, r • • A/ 



(87) 

(38) 

(89) 

(90) 



' S S |-U 

where Uj , iq , etc. denote functions evaluated at the i tn nodal point 

i x. 
i 

on the immersed surface, S 0 . The complex matrices occurring in Fqs. 

Co 7 -90 ) are given by 




J 



7 ; >>) 



'A S 



(91 ) 
(9?) 



( 93 ) 




J 



In evaluation of the integrals in Eqs.(83) and (91-98) it is generally 
adequate to simply evaluate G and its derivatives at the centroid of the 
panel and multiply by its area. However, when point i is either equal to 
j or near j it is necessary to take more care in integrating the 1/R- 
singularity in G. The procedures given by Garrison (1978) were used to 
evaluate the integrals on quadral ateral panels. 

To obtain a solution to the second-order problem it is necessary to 

* 



first evaluate the function f (j,3) required in Eq. (49) and as defined by 
Eq.(50). For this purpose we proceed again numerically by dividing the 
mean free surface (or y = 0 plane) in the vacinity of the body into K area 
elements and evaluate f at the nodal points of these elemental areas. The 
numerical procedure is based on the use of Eq.s(87-94) with y^ = 0 for 
purposes of determination of u^ and its derivatives at points on the free 
surface. The function f is then evaluated at the nodal points by use of 
Eq . (50) . 

The solution for the second-order source strength is obtained numerically 
through the integral equation, Eq. (79). That is, Eq . (79) is written as 

J =!,Z r ..N 

k = t,z r -K 



KiiW £j = h ti i-0( ik (4v> f k 



(95) 



where Qf£ y v) is defined by replacing V by in Eq . (S3) and 
k; = t L !7 X . cosh[a«(li+t/i)]]e (96) 



Once the CK -matrices in (95) hove been evaluated it may be solved for the 
source strength distribution on the immersed surface. 
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The hydrodynamic pressure and resulting forces carried to the second 
order as given in Eq.(59) and (69), respectively, require only the second- 
order potential rather than its derivatives. Once f 2 is known, U 2 may be 
obtained through Eq.(78) which, when written in indicia! form, becomes 










(97) 



where ftijCW) 1S defined by Eq.(91). 

NUMERICAL RESULTS 

A computer code has been developed for the special case of a vertical 

circular cylinder (pile). MacCamy and Fuch's solution was used for the 

★ 

first-order solution and for purposes of evaluating f . The second-order 
solution was calculated on the basis of the distributed singularities as 
given by Eq. (78) . 

The particular example calculation was carried out for a cylinder placed 
in water one radii deep. The dimensionless frequency parameter was varied 
over the range 1.2. The results for the first-order horizontal force 

coefficient, the second-order horizontal force coefficient, steady-state 
force coefficient and phase angles is shown in Figure 3. The results appear 
to indicate that the second-order effects are greatest at the lower freq- 
uency range, i.e., for v*<0.6. 

In Figure 4 and 5 the dimensionless horizontal force is plotted for a 
complete wave cycle for two different frequencies (wave lengths). These 
results which correspond to rather steep wave show a rather sizeable non- 
linear effect. Also, the effect of including second-order effects is to 
shift the maximum force towards the phase of the wave crest. 
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FIGURE 5 HORIZONTAL HAVE FORCE ON A VERTICAL CYLINDER 
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